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Abstract 

Recent high-precision measurements of atomic parity-nonconserving transi- 
tions between the 6S and 7S states of cesium allow for a determination of the 
weak nuclear charge with a precision of 1.3%, providing an improved test of 
the standard model at low energy. Implications for new physics are examined 
in terms of low energy effects on the weak charge, in particular contact inter- 
actions and scalar leptoquark limits. Prospects for further improvements are 
described. 
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I. INTRODUCTION 



Nuclei and electrons are bound by electromagnetic interactions that do not violate parity. 
However the exchange of a Z gauge boson between the nucleus and the atomic electrons is 
parity violating. The dominant contribution arises from the vector coupling to the nucleus 
with the axial-vector coupling to the electrons, which allows normally forbidden transitions. 
The nuclear vector current is conserved and the nucleus acts as a source of the weak charge 
Qw-i which is a linear combination of the Z vector coupling to up and down quarks: 

Q w = -2[C lu (2Z + N) + C ld (Z + 2N)] (1) 

where N is the number of neutrons, Z the number of protons and 

14 12 
Ciu = -- + - sin 2 9 W , C ld = - - - sin 2 W , (2) 

at tree-level. The tree-level formula for the weak charge is modified by radiative corrections 
and by physics beyond the standard model, thus a precise determination of Qw provides 
a test of these corrections and bounds on new physics, complementing the results of high 
energy colliders. Theoretically the structure of cesium is the most accurately known (1%) 
among heavy atoms, as it is an alkali-metal atom and it can be described as a valence 
electron and a closed-shell tightly bound core which is relatively unpolarizable. 

Recently a factor of 7 improvement in the measurement of parity-nonconserving transi- 
tions between the 6S and 7S states of 133 Cs was reported || with the use of a spin polarized 
atomic beam. The nucleus has spin 1 = 7/2 and the total angular momentum of the atomic 
S states is then 3 or 4. The experiment measures both 6S(F = 3) — > 7S(F = 4) and 
65 (F = 4) — >• 7S(F = 3) transitions. The two measured transitions differ because of 
nuclear spin-dependent effects. A linear combination of the two values eliminates nuclear 
spin-dependent contributions, allowing a precise determination of the weak charge. The 
result of is: 

Q w = -72.11 ±0.27 ±0.89 (3) 
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where the first error is experimental and the second one is the consequence of atomic theory 
uncertainty. Contrary to previous data the theoretical uncertainty dominates the error. 

The other linear combination allows a determination of the spin-dependent effects, mainly 
due to the anapole moment |J. However it involves strong interaction uncertainties and I 
will not discuss it in the following. 

Assuming the validity of the standard model one can extract the Weinberg angle from 
Qw, in order to have an idea of the constraints implied by the new cesium data. Using the 
formulae of [|]|§ and the updated analysis of [|7| concerning the hadronic vacuum polarization 



corrections to mixing, in the MS scheme at the scale mz, and with a Higgs mass of 
300 GeV, I find sin 0w( m z)~M~s — 0.2267 i 0.0040. In order to compare this result with 
LEP data one should calculate instead sin 2 # e //, which is the usually quoted value, but the 
difference is tiny || , and one can safely neglect it in the following considerations. From the 
previous determination one can see that in sin 2 @w( m z)-Ms ^ ne error is reduced by a factor of 



two with respect to previous determinations from cesium [10]. Now the error is dominated 
by the atomic physics theoretical accuracy which may improve in the near future to 
few parts in 10 3 . The present error in the siia 2 9\y determination is an order of magnitude 
larger than the one quoted by LEP, sin 2 9 e ff = 0.23200(27) |TT|, but cesium is likely to play 
a more significant role in tests of the standard model in the future. For example assuming 
the atomic theory error to go down to 0.2% and the experimental error to stay the same, 
the induced error in the determination of sin 2 6 W would be 0.0013. 

II. CONTACT INTERACTIONS 

Recently much interest was devoted to contact interactions and scalar leptoquarks as 
they might account for the excess of high Q 2 events at HERA El in e + p collisions. Atomic 



parity violation puts severe bounds on quark-lepton four-fermion contact interactions fll"3"|.|14' 
The relevant lagrangian is written in the form |p~5|: 
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£ 4/= E TT^h^l^i^^ ( 4 ) 

i,j=L,R;q=u,d VHj) 

where €ij = ±1. The contact interaction produces a shift of the weak charge: 

/ d U £Z U £1 U C U \ 

T 2\ I t RL t LL e LR, , e RR \ 



6Q W = — (143TeV 



\u 2 An 2 \u 2 \u 2 
s IV RL 1V LL IV LR ly RR , 



- (160TeV 2 ) ( - - "H + -^2 1 (5) 

and using the experimental result for the weak charge 0] and the theoretical standard model 
value —72.88 ± 0.06 of [|10J, eq. @ gives the limits of table I, which are larger than those 



coming from high-energy colliders (see for example |16| and recent bounds by CDF |17|| and 



LEP 2 One should be careful, when comparing table I with the bounds reported in the 

literature, that the definition of Af^ may differ by a factor of \f^n. The limit coming from eq. 
(||) can however be eluded if the contact interactions are parity conserving or cancellations 
occur (see the discussion in HTDl). 



III. SCALAR LEPTOQUARKS 

Let us now consider the low energy effects of scalar leptoquarks that couple chirally and 
diagonally to the first generation. "Chiral" means that the leptoquark is coupled either 
to left-handed or to right-handed quarks but not to both, in order to avoid unacceptable 
deviations from lepton universality, for example in 7r — > eu, while "diagonal" means that the 
leptoquark couples only to a single leptonic and quark generation (at least approximately as 
the CKM matrix induces mixing for left-coupled leptoquarks, and as a consequence flavour 
changing neutral currents). 



I follow here the notation of [50], where previous cesium data was used to set bounds 
on the leptoquarks. From the general effective lagrangian satisfying baryon and lepton 
number conservation, with the most general dimensionless and SU(3) C x SU(2)l x U(1)y 
invariant couplings one has scalar leptoquarks with the following quantum numbers with 
respect to SU(3) C x SU(2) L x U(l) Y - 5(3*, 1, 1/3), S(3*, 1, 4/3), T(3*, 3, 1/3), L>(3, 2, 7/6), 



.0(3,2,1/6). S and D couple only to right-handed quarks, T only to left-handed quarks, 
while S and D can couple to both. 

In a class of supersymmetric theories it is natural to have Sl, S r of type S and Dl, D r 
of type D P2"| , where the subscripts i? and L refer to the helicity of the quark coupled to the 
leptoquark. The requirement of having chiral couplings in the sense indicated before may 
therefore be fulfilled, if mixing between the left- and right-coupled leptoquarks is small. 

In the following I shall be mainly interested in Dl, D r and D leptoquarks, which are 
relevant in the discussion of the excess of high Q 2 events at HERA in e + p collisions. Their 
Yukawa interactions are given by 

£ = g R (eu R D R ' 5/3) + 9u R D^) + g & (ed R D^ + ud R D^) 

+ g L (eu L D L - Fj/3) + ed L D { L - 2/3) ) + c.c. (6) 

where for the Dl I neglected second and third generation couplings, as the leptoquark is 
assumed to be diagonal to a good approximation. The superscript in parentheses in the 
previous formula refers to the leptoquark electromagnetic charge. 

One can easily derive the effective low energy 4-fermion interaction due to leptoquark 
exchange and obtain the corresponding shift in the weak charge. In order to have a compact 
notation, it is useful to introduce the parameter rji which has the value 1 when the leptoquark 
I is present and otherwise is zero. The contribution to Qw is: 

5Q ™ = ~272G~ F (m) 2 + VSr ~ VDl + VDr ~ VT)i2Z + N) 

+ (Vs-Vd l +Vd-2Vt)(Z + 2N)} (7) 

where Gf is the Fermi constant. Using the experimental result for the weak charge 0] and 
the theoretical standard model value —72.88 ± 0.06 of |TIJ, one obtains the limits of table 
II, in terms of the ratio leptoquark mass M over coupling g. 

Limits from other experimental data can further constrain the leptoquarks, depending 
on the nature of the leptoquark, for example flavour changing neutral current limits or 
universality in leptonic it decay. Moreover note that the limits coming from eq.(^) are 
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valid when one assumes that there is only one leptoquark multiplet and that there is no 



mass splitting within the multiplet (for a detailed discussion see appendix B of P0|). The 
presence of more than one leptoquark multiplet can both improve or reduce the bound, as 
different leptoquarks may contribute to the weak charge with the same or opposite sign. 

I will not analyse the implications of leptoquarks for the high-Q 2 HERA data in detail, as 
this is the subject of dedicated papers ^3|. Note however that Dl,r and D leptoquarks are 



produced in e + p collisions in the s-channel, whereas they contribute to e~p non-resonant ly. 
If a mass of 200 GeV is assumed for the leptoquarks, a bound on the couplings can be 
derived, 

g L , R < 0.09 g b < 0.08 (8) 

at 95 % C.L. while the coupling required to explain the HERA anomaly is lower, of the 
order of 0.02 for a scalar leptoquark coupled to u quarks (like Dr), and 0.04 if coupled to d 
quarks (like D) and thus within the permitted region. 

Limits form cesium can be combined with direct leptoquark searches to perform a two 
parameter fit in the mass and coupling of the leptoquarks. In the following \ 2 analyses a 
small number of events has to be compared with theory. I use for them the method of least 
squares with a \ 2 function for Poisson-distributed data, which asymptotically behaves like 
a classical x 2 [HI- The goodness-of-fit (confidence level) is evaluated approximatively as if 
data were Gaussian-distributed. 

Leptoquark production at Tevatron is insensitive to the g coupling (the only requirement 
is g > 10~ 12 due to the event reconstruction algorithms used in the experimental analysis), 
but detection of the leptoquark decay products depends on the leptoquark branching fraction 
to £q and vq. A recent analysis from the DO collaboration was used in the fit |^5| . The leading 
order parton production cross-section [^6j was convoluted with parton distribution functions 
27], at a scale \i = M (this choice of the scale at leading order is supported by the next- 



to- leading result of |28||), in order to produce the theoretical cross section at Tevatron. The 
result of the fit is shown in Fig.l, as a 95% C.L. bound on the mass and coupling of the D 



scalar leptoquark (limits for Dl and Dr are similar and are not shown). It is assumed that 
D has 100% branching fraction to eg. Fig. 2 shows the same limits assuming 50% branching 
to eq and 50% to v e q. 

Assuming that the excess of events seen at high-Q 2 in e + p collisions at HERA |12j is due 
to the production of a scalar leptoquark, I calculated bounds in the plane (M,g) using cesium 
data || and the leading order parton cross-section for scalar leptoquark production [21] with 



parton densities [27], taking into account initial state photon radiation from the positron. 



The combined integrated luminosity of 34.3 pb _1 from HI and ZEUS collaborations at \/s = 
300 GeV in the e + p — > eX mode was used, taking into account that the two experiments 
have seen 24 events with Q 2 > 15000 GeV 2 , while the standard model expectation is 13.4 ± 1 
events. Fig. 3 refers to a D leptoquark, while Fig.4 to a Dr leptoquark. The limits for a Dl 
leptoquark assuming that the members of the doublet are degenerate in mass, are similar to 
those of Fig.4 and are not shown. No information was given to the fit on the invariant mass 
distribution of the anomalous HERA events, in order to see if the preferred (M,g) values 
were consistent with those required by the leptoquark interpretation of HERA anomalous 
events. They turn out to be consistent, however x 2 is small in a narrow band over a wide 
range of masses. Note that if the contour plot is performed using a too coarse-grained grid 
of points, this can be missed. One can at most say that for the D leptoquark the preferred 
region is for M < 220 GeV, g < 0.06 (for a mass of 200 GeV the preferred value of g is 
0.034) and for the D R leptoquark M < 250 GeV, g < 0.07 (for a mass of 200 GeV the 
preferred value of g is 0.017). The preferred region for the fit is obtained demanding x 2 < 2, 
note however that for low x 2 the Gaussian confidence-level estimate is not appropriate. 
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TABLES 





e u = 1 


e u = -1 


e d = 1 


e rf = -1 


LL, LR 


7.4 


11.7 


7.9 


12.3 


RR, RL 


11.7 


7.4 


12.3 


7.9 



TABLE I. 95% C.L. lower bounds from contact interactions on A?- in TeV units assuming o 
one of the ef • terms is responsible for the shift in the weak charge. 





S L 


S R 


D L 


D R 


T 


S 


D 


M/g 


1500 


2300 


2200 


2300 


2700 


2500 


2500 



TABLE II. 95% C.L. lower bounds on scalar leptoquarks in GeV units for the ratio M/g. 



10 



FIGURES 
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Fig. 1 - 95% CL. limit in the mass M (GeV) and coupling g of a D leptoquark from 
atomic parity violation in cesium || and the direct search at TEVATRON (data from the 



DO collaboration [25]), assuming a branching ratio of the leptoquark to ed equal to 100%. 
The allowed region in the parameter space is the grey area. Darker regions correspond to 
lower x 2 values. 
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Fig. 2 - 95% CL. limit in the mass M (GeV) and coupling g of a D leptoquark assuming 
a branching ratio of the leptoquark to ed of 50% and to v e d of 50%. 
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Fig. 3 - 95% C.L. limit in the mass M (GeV) and coupling g of a D leptoquark from the 
combination of cesium data and the cross-section excess at HERA. 



0.2 



0.15 



. 1 



0.05 








140 160 180 200 220 240 



Fig. 4 - 95% C.L. limit in the mass M (GeV) and coupling g of a D R leptoquark from the 
combination of cesium data and the cross-section excess at HERA. 
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